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Abstract. We present a brief summary of the main results from our multi- 
dimensional, time-dependent simulations of gas dynamics in AGN. We focus on 
two types of outflows powered by radiation emitted from the AGN: disk winds 
and winds driven from large-scale inflows. We show spectra predicted by the 
simulations and discuss their relevance to observations of broad- and narrow-line 
regions of the AGN. We finish with a few remarks on whether these outflows can 
have a significant impact on their environment and host galaxy. 



1. Introduction 

While AGN are very strong sources of electromagnetic radiation, they are also 
sources of outflows of matter and magnetic energy. There is a growing body of 
evidence that outflow s are quite common and are an integral part of AGN (e.g., 
Crenshaw et al. 20021 : Richards Hall! [2004). Moreover, radiation and mass 



outflows in luminous sources are physically coupled as matter can absorber, 
scatter and emit radiation and radiation can affect dynamics and ionization of 
matter. For example, powerful mass outflows in quasars are very likely winds 



matter, tor example, powertul mass outflows m quasars are very likely winds 
driven by radiation from accretion disks (e.g., reviews by IKonigll 120061 : IProgal 



2007al h In addition, as we will illustrate here, the same radiation can drive 



outflows from large scale inflows. 

In this paper, we briefly summarize results from numerical simulations of 
disk winds and large-scale outflows. We review the properties of these outflows 
and compare them with those observed in AGN. This will allow us to assess the 
level of our understanding of AGN outflows and their potential role in the AGN 
feedback. 



Radiation-Driven Outflows from Accretion Disks 



The successes of modelling outflows driven by radiation pressure on sp ectra 



lines (line driving) from OB stars (e.g.. ICastor et al.lll975 



Puis et al 



Progal 120051 ) motivate 



2008) and 



from accretion disks in cataclysmic variables (CVs; e.g.. 
applica tions of a similar p hysics to model outflows i n AG N. 

In iProga et al. ( 20001') and Proga fc Kallmanl ( 20041 ). we adopted the ap- 
proach from Proga et al. ( 19981 . 19991 ) to calculate axisvmmetric time-dependent 
hydrodynamical models of line-driven winds from accretion disks in AGN. To 
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Figure 1. A map of the poloidal velocity field of the radiation-driven disk 
wind model described in the text. The rotation axis of the disk is along the 
left hand vertical frame, while the midplane of the d i sk is along the lower 
horizontal frame. The figure is from lProga fc Rahman! (|2004l) . 



apply the line-driven disk wind model developed for CVs to AGN, we took into 
account, in an approximated way, the difference in the spectral energy distribu- 
tion between AGN and CVs. In particular, we introduced three major modifi- 
cations to our previous approach: 1) calculation of the parameters of the line 
force based on the wind properties, 2) effects of optical depth on the continuum 
photons, and 3) radiative heating and cooling of the gas. 

Our AGN wind calculations follow (i) a hot and low density flow with 
negative radial velocity in the polar region, (ii) a dense, warm and fast equatorial 
outflow from the disk, and (hi) a transitional zone in which the disk outflow is 
hot and struggles to escape the system. Fig. 1 shows an example of a velocity 
field of a AGN disk wind model. 

To make a direct comparison with observations, w e calculated synthetic lin e 
profiles based on our model for the C IV A 1549 A line ( Proga &; Kallman 20041 ). 
The synthetic line profiles show a strong dependence on inclination angle: the 
absorption forms only when an o bserver looks at the CE t hrough the fast wind 
(i.e., i > 60° as shown in Fig. 3 of lProga Sz Kallmanl 12004 : see also Fig. 2 here). 
This i dependence can explain why only 10% of QSO have BALs. 

The model also predicts high c olumn densitie s for t he inclination angles at 
which strong absorption lines form. Schurch et al. ( 20091 ) co mputed broad band 
spectr a for various inclination angles using the simulation from lProga fc Kallman 
(12004T ). ?ig. 3 presents some examples and shows that the model is consistent 
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Figure 2. Theoretica l profiles of Civ i549A based on the wind model of 
iProga fc Kallmanl (|2004l ) as a function of inclination angle, i (see bottom right 
corner of each panel for the value of i). The solid lines show the profiles due 
to the absorption only whereas the dashed lines show the profiles due to the 
absorption and emission (the line source function includes resonance scatter- 
ing and thermal emission). Note the sensitivity of the lines on the inclination 
angle. The zero velocity corresponding to the line center is indicated by the 
vertical line. 



with the observations, i.e., BAL QSO a re underluminous in X-rays compared t o 
their non-BAL QSO counterparts (e.g-. lBrandt et al1l2000l : iGiustini et al.ll2008l ). 



3. Radiation-Driven Outflows from Large-Scale Inflows 

In the previous section, we showed that the powerful radiation from an AGN 
can drive a strong wind from an accretion disk, the place where the radiation 
is produced. In turn, this wind has a significant impact on the radiation. How- 
ever, the same radiation can also change the dynamics of the material at large 
distances from the radiation source. 

We h ave begun studying gas dynamics in AGNs on subparsec and pa r- 
sec scales dProgal l2007bl : iProea et all 120081 : iKurosawa Proeal 120081 . l2009al lbl: 
Kurosawa et al. 20091 ). In the following, we briefly summarize the main results 
of this work. 

In iProeal (|2007bl ). we calculated a series of models for non-rotating flows 
that are under the influence of super massive BH gravity and radiation from an 
accretion dis k surrounding the B H. Generally, we used the numerical methods 
developed by IProga et al. ( 2000l ). Our numerical approach allows for the self- 
consistent determination of whether the flow is gravitationally captured by the 
BH or driven away by thermal expansion or radiation pressure. 

For a 10 8 M BH with an accretion luminosity of 0.6 of -^Edd (the same 
parameters as in the disk wind simulations), we found that a non-rotating flow 
settles quickly into a steady state and has two components (1) an equatorial 
inflow and (2) a bipolar inflow/outflow with the outflow leaving the system 
along the pole. The first component is a realization of Bondi-like accretion flow. 
The second component is an example of a non-radial accretion flow becoming an 
outflow once it is pushed close to the rotational axis of the disk where thermal 
expansion and radiation pressure can accelerate the flow outward. The main 
result of these simplified calculations is that the existence of the above two 
flow components is robust yet their properties are sensitive to the geometry, 
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Figure 3. The 0.3-13 keV transmitted X-ray spectra (E Le) from five dif- 
ferent lines-of-sight, with different inclinati ons, through the hydrody namic 
simulation of an AGN outflow presented in iProga fc Kallmanl (|2004j ) . The 
spectra were calcu lated with XSCORT (v5.18) using snapshot 800 of the 
IProga fc Kallmanl simulation to provide self-consistent physical outflow prop- 
erties. The input ionizing power-law spectrum is shown in black. The spectra 
correspond to inclinations of 8 = 50, 57, 62, 65, & 67°, (top to bottom). 
These l.o.s were chosen to highlight the range of spectral shapes that result 
from the diff erent physical proper ties throughout the simulated wind. The 
figure is from lSchurch et all (l2009h . 



SED of the radiation field, and the outer boundaries. In particular, the outflow 
power and the degree of collimation are higher for the model with radiation 
dominated by UV/disk emission than for the model with radiation dominated 
by X-ray/central engine emission. This sensitivity is related to the fact that 
thermal expansion drives a weaker and wider outflow, compared to the radiation 
pressure. 

Rotation of the inflowing gas changes the geometry of the flow because the 
centrifugal force p revents gas from reaching the rotational axis (see Fig. 1 in 
Proga et al.1 12008). This, in turn, reduces the mass outflow rate because less 



gas is pushed toward the polar region. We also found that rotation can lead to 
fragmentation and time variability of th e outflow. As the flow fragments, cold 
and dense clouds form (Figs. 4 and 12 in Proga et al. 20081) . 



Three-dimensional effects are also important. In Kurosawa Prog a (2008), 
we considered effects of radiation due to a processing accretion disk on a spherica l 
cloud of gas around the disk. On the othe r hand , in Kurosawa Sz Proga ( 2009ah . 
we recalculated some models from papers iProeal (|2007bT ) and lProga et al.l "(j2008) 
in full 3-D. Our 3-D simulations of a nonrotating gas show small yet noticeable 
nonaxisymmetric small-scale features inside the outflow. However, the outflow 
as a whole and the inflow do not seem to suffer from any large-scale instability. 
In the rotating case, the nonaxisymmetric features are very prominent, especially 



Figure 4. Three-dimensional hydrodynamical simulations of outflow forma- 
tion via redirection of accreting gas by the strong radiation from an accretion 
disk around a super massive black hole with its mass Mbh = 10 8 Mq. The in- 
falling gas is weakly rotating (sub-Keplerian) , and the Eddington ratio of the 
system is 0.6. The volume rendering representation of the density distribu- 
tions is shown. The outflow morphology is bi-conical, but the flow contains rel- 
atively cold and dense cloud-like structures whic h resembles those ob s erved i n 
the NRLs of Seyfert galaxies. The figure is from Kurosawa fc Progal ([2009b). 



in the outflow which consists of many cold dense clouds entrained in a smoother 
hot flow (e.g., see Figs. 4 and 5). The 3-D outflow is nonaxisymmetric due to the 
shear and thermal instabilities. Effects of gas rotations are similar in 2-D and 
3-D. In particular, gas rotation increases the outflow thermal energy flux, but 
reduces the outflow mass and kinetic energy fluxes. In addition, rotation leads 
to time variability and fragmentation of the outflow in the radial and latitudinal 
directions. The collimation of the outflow is reduced in the models with gas 
rotation. The main difference between the models with rotation in 3-D and 2-D 
is that the time variability in the mass and energy fluxes is reduced in the 3-D 
case because of the outflow fragmentation in the azimuthal direction. 

To be able to compare these new simulations with observations, we are in a 
process of computing synthetic line profiles, broad band spectra and maps. How- 
ever, even without these diagnostics we can check if the models are consistent 
with the data. F or example, we can compare the kinematics study of NGC 4151 



Das et ahl ( 20051 ) with the velocity of the cold clouds (cf. Fig. 5) projected (v 



pro) J 



toward an o bserver at the in clination angle i = 45° , which is the inclination of 
NGC 4151. Das et ahl ( 20051 ) used the kinematics model of the outflows with a 



bi-conic radial velocity law, and found a good fit to their observations when the 
opening angle of the cone is ~ 33°. Interestingly, we find the opening angle of 
the outflows in our is also about 30° (cf. Figs. 5). 
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Figure 5. Spatial distributions of the "cold clouds" in the model shown 
in Fig. 4. The grayscale image shows the density map of the cold clouds in 
logarithmic scale (in cgs unit) on the z—x plane. The cold clouds here are 
defined as the gas with its density higher than p m - ln — 1.6 x 10~ 20 gcm~ 3 and 
its temperature less than T max = 1.6x 10 5 K. The clouds are not spherically 
distributed, but located near the bi-conic surface (which appears as an X- 
shaped pattern here) defined by the outf lowing gas. Note that the length 
scale are in units of pc. The figure is from Kurosawa fc Progal (|2009bl ) . 




Figure 6 shows i> pro j of the clouds plotted as a function of the projected 
vertical distance, which is the distance along the z-axis in Fig. 5 projected onto 
the plane of the sky for an observer viewing the system with i = 45°. The 
figure shows that the clouds are accelerated up to 250 km s^ 1 until the projected 
distance reaches ~ 4 pc, but the velocity curve starts to flatten beyond this 
point. Towards the outer edges (near the outer boundaries), the curve begins to 
show a sign of deceleration, but not so clearly. We note that the hot outflowing 
gas, on the o ther hand, does show dece leration at the larger radii in our models 
(cf. Fig. 9 in Kurosawa Sz Progal 12009a). Al t houg h the physical size of the long 
slit observation of NGC 4151 bv iDas et all (|2005T l is in much lager scale (~ 50 
times larger) than that of our model, their radial velocities as a function of the 
position along the slit (see their Figs. 5 and 6) show a similar pattern as in our 



model (Fig. 6). The range of t> pro j in our model is about 



while the range of the observed radial velocities in IDas et al 
—800 to 800 km s -1 , which is comparable to ours. 



250 to 300 km s" 1 
(|2005l ) is about 



4. Concluding Remarks 

The models presented here, which numerically simulate the outflows driven by 
radiation from AGN, are in many respects consistent with observations. Clearly 
more work is needed to test the models and improve them. However, the first 
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Figure 6. The velocities of the cold cloud elements (as in Fig. 5) projected 
toward an observer located at the inclination angle i = 45° are shown as 
a function of the projected vertical distance (the distance along the z-axis 
in Fig. 5, but projected on to the plane of the sky for the observer viewing 
the system with i — 45°). The negative projected distance indicates the 
clouds are found in the low er half of the projection plane. The figure is from 
Kurosawa k. Progal (|2009bh . 



few steps toward the development a self-consistent physical model of the AGN 
outflows has been taken. 

Applications of the model are not limited to AGN because other astro- 
physical objects - such as X-ray binaries, in particular micro-quasars - have a 
similar geometry and can be understood within a similar physical framework. 
In addition, having a physical model of the AGN outflows we can apply it to 
the so-called AGN feedback problem. 

Some results from the outflow models have been al r eady incorporated to 
galaxy evolution calculations. In particular, ICiotti et al.l ( 20091 ). improved and 
exten ded the accretion and feedback phys ics explored in their previous papers 
(e.g., ICiotti k Qstrikerl Il997l . l200ll. 120071). By usin g a high-resolution one- 
dimensional hydrodynamical code, Ciotti et a l. (2009) studied, the evolution of 
an isolated elliptical galaxy, where the cooling and heating functions include pho- 
toionization and Compton effects, and restricting to models which include only 
radiative or only mechanical feedb ack (the latter, in the form of disk winds pro p- 
erties of which were adopted from Proga et al. 2000] ; Proga Kallman 20041 ) . 

These recent calculations confirmed that for Eddington ratios above 0.01, 
both accretion and radiative outputs are forced by feedback effects to be in burst 
mode, so that strong intermittencies are expected at early times, while at low 
redshift the explored models are characterized by smooth, very sub-Eddington 
mass accretion rates punctuated by rare outbursts. However, the explored mod- 
els always fail some observational tests. For the high mechanical efficiency of 
10 -2,3 as adopted by some investigators, it was found that most of the gas 
is ejected from the galaxy, the resulting X-ray luminosity is less than a value 
typically observed and little super massive black hole growth occurs. However, 
models with low enough mechanical efficiency to accommodate satisfactory black 



8 



hole growth tend to allow too strong cooling flows and leave galaxies at z = 
with E+A spectra more frequently than is observed. Surprisingly, it was also 
found that both types of feedback are required. Radiative heating over the inner 
few kilo parsecs is needed to prevent calamitous cooling flows, and mechanical 
feedback from active galactic nucleus winds, which affects primarily the inner 
few hundred parsecs, is needed to moderate the luminosity and growth of the 
central black hole. Models with combined feedback pass more observational 
tests (Ciotti et al., in preparation). Thus it emerges from these simulations that 
to solve the AGN feedback problem/explain all aspects of the galaxy evolution, 
more than one form of feedback is needed. 
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